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Abstract - We address an unresolved issue in the physics of low-dimensional many-body systems: the ques- 
tion of whether or not a random field can produce order at low temperatures for statistical mechanical systems 
possessing continuous internal symmetries. Concretely, we verify that the XY model in a uniaxial random 
field orders in two and three dimensions. The direction the system orders is peipendicular to the randomness 
for any choice of symmetry breaking field with nonzero projection peipendicular to the randomness. The 
result is particularly relevant in two dimensions, where there are a number of competing effects - quasi-long 
range order of the pure system and strong fluctuations of the random field. While we consider only classical 
systems explicitly, the effect is robust and our work has implications for quantum systems as well, producing 
ordered phases any dimension. 



Introduction. It is generally understood that disorder plays an 
important role in the behavior of one-particle and many-body 
physics. The most well-known example stems from the work 
of Anderson |J6l on transport properties of electrons in crystals 
with impurities; disorder can localize electron wave functions 
making conductivity negligible. Other important examples of 
the effects of disorder include the "rounding" of first order tran- 
sitions in low dimensional classical equilibrium systems 15.24.1 . 
This work was recently extended to quantum systems in ||4]. 

The rounding effect aside, the behavior of interacting many- 
body systems in the presence of disorder is less well under- 
stood. For example, it was only recently demonstrated by an in- 
depth analysis that Anderson localization persists in the pres- 
ence of weak interactions between electrons dTJ. Localization 
of Cooper pairs was recently observed in the vicinity of the 
BCS superconducting transition due to the presence of disorder 



In these examples, as well as in subjects such as quantum 
computing, the existence of quenched randomness in a system 
is viewed as an unwanted effect which disrupts ordering and 
coherence. It is less well known that randomness can itself 
create ordering for systems in which the order parameter has a 
continuous symmetry. We refer to this effect below as random 
field induced order (RFIO). 

The basic model for RFIO is an XY model in a uniaxial ran- 
dom field. It is a classical equilibrium statistical mechanics 
model in which spins take values on the unit circle and interact 
in a ferromagnetic way. There is additionally a random field 
acting in the vertical direction only. In more precise terms and 



generalizing the setting, we will below consider the following 
family of models. Let Ajv = {-A^/2, ... ,iV/2 - 1}''. For 
X E An, spin variables ax lie in the unit sphere §"^^ C M", 
n > 2. We define the (random) Hamiltonian via 
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The first sum is over adjacent pairs x,y in An while the final 
sum consists of vertices at the boundary of A^v. The ax's form 
a family of independent, identically distributed /c-dimensional 
standard Gaussian vectors, where we have taken k < n. The 
last term is for symmetry breaking purposes. When defining 
Gibbs measures and states, the a priori measure on spin space 
is uniform surface measure on §"~^. We refer below to these 
models collectively as RFO(n; k) models. In this notation, the 
XY model with a uniaxial field is the RFO(2; 1) model. The 
latter will be shortened to the RFO(2) model. 

The RFO(2) model was first studied in the 1970's via mean 
field theory and a renormalization group calculation [|2]. The 
system appears to order, but only perpendicular to the axis 
on which the randomness acts. In low dimensions, particu- 
larly on the two-dimensional square lattice Z^, the mean field 
and renormalization group calculations are less reliable, as the 
conclusions run counter to other effects, such as the Mermin- 
Wagner theorem and the Imry-Ma argument that strong fluctu- 
ations disrupt ordering ll24l (more about the relation between 
RFIO and these results appears below). Thus in the 1980's 
multiple groups addressed two-dimensional behavior by more 
specialized renormalization group analyses llT2][T3]|26l . One 
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group lfT2l concluded there is a low temperature paramagnetic 
phase while the other ll26l concluded there is an intermediate- 
temperature ordered phase from which they extrapolate the 
low-temperature behavior Interesting, but only tangentially re- 
lated, rigorous work was done in the 1990's on ground states in 
the strong field regime in II16II17I . 

Both analyses of the RFO(2) model in a weak field are prob- 
lematic. The work 01211131 is based on an approximation via a 
random field sine-Gordon model, which is extremely unstable 
with respect to the noise (see our related discussion of the ran- 
dom field Gaussian model below). The paper ||26l suggests the 
ordered phase persists at a fixed temperature even as the field 
strength e is taken to 0. Such a conclusion is dubious as sim- 
pler versions of the model, such as a model which replaces the 
random field by a periodic field with period two has a critical 
inverse temperature which diverges as the field strength tends 
to in two dimensions. 

With this background it seems useful to provide a definitive 
analysis of the the RFO(2) model in two dimensions. We do 
this and more in two main results below, the method of ver- 
ification being a rigorous proof. In Theorem [T] we state, for 
dimension d ~ 2,3, that ordering in the direction perpendicu- 
lar to randomness does occur at low temperature. Moreover we 
show that this is the only direction in which ordering occurs in a 
strong sense: it is independent of the direction of infinitesimal 
field used to break symmetry. This is part of a more general 
result for RFO(n; n — 1) models. A second result. Theorem 
121 deals with RFO(n; k) models for any k < n and in any 
dimension d > 2: The projection of the spin variable at any 
fixed vertex onto the subspace which supports the randomness 
is necessarily small. 

Let us remark that the issue of RFIO resurfaced over the last 
five years in the work of three separate groups. With the ad- 
vent of controllable interacting Bose-Einstein condensates in 
optical traps, the effect was suggested as a possible response to 
the presence of certain kinds of experimentally realizable dis- 
order II30II31I . Here the (pseudo-)spin variables arise from in- 
ternal structure of the atoms, the tuning of interactions and the 
structure of the optical lattice. Building on 123], the same sort 
of mechanism was put forward as the reason for the splitting 
of Landau level degeneracy in experiments on graphene fl], 
though the particulars of this model, namely the fact that the 
derived random field is the curl of random coupling constants, 
means that the fluctuations of the disorder are subcritical from 
the perspective of the Imry-Ma argument. Finally, in [15], this 
question was posed during investigations of whether real-space 
coarse-graining procedures preserved a version of the spatial 
Markov property which characterizes Gibbs measures. 

To explain the prior controversy and subtlety of our results 
in two dimensions we recall, for the convenience of the reader, 
the behavior of related models. 

The Pure 0(2) Model. In our language this is the case we take 
n = 2 and e = 0, so there is no randomness. When d = 2 
all Gibbs states are rotationally invariant in the thermodynamic 
limit (this is the content of the Mermin- Wagner theorem, see 
ifTTllZTl among many other works) and, in particular, there is 
no residual magnetic order There is however a Kosterlitz- 



Thouless phase transition B 191125 1 expressed by a change in the 
behavior of the spin-spin correlation function {(Jx ■ cfy) from 
exponential to algebraic decay in terms of jcc — y|. If d > 3 
residual magnetic ordering does occur 11201 . 
The Random Field Ising Model. In this case we constrain the 
spins to point only in the vertical direction, the same as the 
random field. When d > 3 and weak disorder, it turns out 
that residual magnetic ordering persists |l9]l22]. In dimension 
d = 2 there is no magnetic order for any inverse temperature 
/? and any strength e, [Sj. This is the rounding effect alluded 
to in the introduction. Physically, see ll24l . the reason for this 
is that at all scales there are random fluctuations strong enough 
to overcome the loss in energy due to mismatched spins on the 
interior and exterior of domains. 

The Random Field Gaussian Model. The last model we wish 
to mention is obtained by replacing the vector valued spins 
(7x G §^ by a field (j)^ € M. and otherwise retaining the setup 
described above. This model appears, among other places, in 
ll8l lT4l . In a finite volume A with boundary conditions, we let 
— Aa denote the discrete Laplace operator on A with Dirichlet 
boundary conditions. We have 
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Here and below E denotes the average over the randomness in 
the model of interest. 

If e > and if d < 4, the second term on the RHS grows 
with A while when d = 2, it grows even after taking gradients 
in both arguments x,y. 

If one believes in magnetic ordering for the RFO(2) model 
in two dimensions, the last example in particular should give 
pause as this is precisely what is obtained by expanding the 
RFO(2) Hamiltonian ^ in angular coordinates around a fixed 
direction. It suggests that either the ansatz of having order is 
flawed or that the fact that spins are constrained to lie on §^ 
plays a crucial role in keeping the influence of field fluctua- 
tions under control. This effect is presumably the reason for 
the paramagnetic low temperature phase prediction in lll2l . On 
the other hand it seems to have been neglected in the arguments 
presented in II26II31I . 

Order-by-Disorder. The phenomenon of "Order-by-Disorder" 
provides a class of systems which exhibit ordering due to vari- 
ous types of fluctuations. The most relevant example, first con- 
sidered by Henley ll2Tl . concerns a model Hamiltonian on 1? 
of the form 
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with I Ji| < 2J2. The ground-states for this (frustrated) system 
are obtained by choosing a purely anti-aligned configuration of 
spins on each of the even and odd sub-lattices of 1? and are 
thus parameterized by two angles: an angle between the spin 
at (0, 0) and the ei-axis and relative angle between the spin at 
(0, 1) and (0, 0). The degeneracy of ground-states is partially 
lifted under the introduction of two types of "disorder". The 
first type consists in passing from to positive temperature. 
More relevant to the RFIO is a second mechanism: site dilution. 
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Vertices of 1? are deleted from the system independently with 
probability p ^ 1. According to the calculations in lISTl . at 
and low temperature the system prefers the ground states with 
relative angle between spins at (0, 0) and (0, 1) to be fixed at 

The crucial difference between this site diluted model and 
the RFO(n) models we consider concerns fluctuations due to 
randomness. These are substantially weaker in the site diluted 
order by-disorder setting. In particular, there is an analog to the 
field (?x introduced below but the fluctuations of this field are 
about as singular as the four dimensional version of the RFO(2) 
model. While the site diluted order-by-disorder problem has 
not been rigorously addressed, this feature suggests the conclu- 
sions in II2TI are reliable. It would be interesting to see if our 
methods can be adapted to this case. 

Main Results: In the following statements, P denotes the pro- 
jection operator onto the k dimensional subspace supporting 
the distribution of a^- 

Recall ([T]i and let (•)^ " denote the corresponding random 
Gibbs state on A^v at inverse temperature f3. For the statement 
of our main result, let us define the block average observables 
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Theorem 1. Let d € {2,3} be fixed and consider the RFO(n;n- 
1) model on 11^. Let e„ he a fixed unit vector so that P • e„ = 0. 
There is an interval (0, eo) so that the following holds: 

For any e G (0, eo) there are ^(e) > O|loge|^^" and 
l3o{e,0 > Oso that if 



u e S"~\u-e„ > ^ 



and j3 > /?o, we have 
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for any z e U^. Moreover ^(e) — > a.? e — > 0. 

A similar effect was previously observed in XY models with a 
weak uniaxial field alternating direction in a chessboard fashion 
in classical ll2l [T5l and quantum |3j models. 

The condition u ■ e„ > ^ is technical; the result should hold 
for any choice u • e„ > independent of e as long as e is small. 
It is possible that further analysis will remove this condition. 
The proof of Theorem [T] also gives lower bounds on /?o(e) of 
order e~^. The analysis of reference |56l suggests the critical 
temperature is e-independent. Resolving this discrepancy is an 
interesting open issue. 

Even at temperature. Theorem [T] strengthens the picture 
presented in 12611311 . Our methods show that, as e tends to 0, in 
the ground state the length of the projection of most spins onto 
the direction perpendicular to the randomness tends to 1. 

Intuition and |2j suggest similar statements should hold at a 
rigorous level for d > 4 and moreover that the effects of large 
fluctuations should be easier to control. Our technical estimates 
bear this out, but currently we cannot provide a complete proof 



in this case (see the discussion below). This is connected with 
the fact that control of the Dirichlet energy of functions pro- 
vides better point-wise control in low dimension than in high 
dimension. 

Next we present a weaker result for general RFO (n; fc) mod- 
els in any dimension d> 2. Let ctq be a fixed spin configuration 
on Z'' and let {■)'^ denote the Gibbs state on A^v with boundary 
condition a^ 

Theorem 2. Let d > 2 and k < n be fixed and consider the 
RFO{n; k) model on 1/^. There is an interval (0, eo) so that for 
any boundary condition cjq the following holds: 

For any e G (0,eo) there exists a I3q{€) so that for the 
RFO{n; k) model witheand /3 > /3o 
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For d e {2, 3, 4}, Theorem|2]is close in spirit to Theorem 4.4 
in the paper by Aizenman and Wehr J5l|, though neither result 
implies the other 

We make no rigorous statements regarding magnetic order- 
ing in RFO(n; k) models for k < n — 2, although ordering 
should occur if d > 3. If d = 2 we expect the RFO(n; n — 2) 
models have a Kosterlitz-Thouless transition and that if fc < 
n ~ 2, the model exhibits exponential decay of correlations 
for all f3 in correspondence with the expected behavior of pure 
0(71 — fc) models when d = 2. 

The choice of standard Gaussian vectors is not crucial. The 
proof of Theorem [T] applies for variables with sub-Gaussian 
tails in d = 2 and for any distribution having arbitrary expo- 
nential moments if d = 3. The proof of Theorem |2] requires 
even weaker assumptions on the tail behavior 
Extensions: The effect described in this note should be very ro- 
bust, applying to classical and quantum systems and both sym- 
metric and biased disorder Examples include: 

Biased Disorder: Suppose the cJa. = X]"=i "^^l*) ^i'^h 
dx{i) independent Gaussian vectors with non-zero means 
E[aa;(i)] = I'i. For example, if the Vi from an orthonormal 
basis for R" ordering will occur along the directions J^i i^i- 
Such a model should be treatable by our methods combined 
with Pirogov-Sinai theory. 

Quantum models : Whenever the quantum system without 
disorder may be expressed in terms of pseudo-spin order pa- 
rameters with continuous symmetry, effects analogous to those 
described above are possible. Examples include Bose and 
Fermi Hubbard models as discussed in ll30l and graphene with 
fluctuating interaction strengths IT] . 

Proofs. For the interested reader, we now give a brief sketch 
of the proof. For simplicity, we restrict the discussion to the 
RFO(2) model. Demonstration of the remainder of Theorem 
[1] follows along the same lines. The proof of Theorem |2] is 
actually simpler because we do not need to distinguish different 
phases. For notational purposes set 
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We use the notation \R\ for the cardinahty of a finite set R C 

Calculations for Small Boxes. To understand the issues we 
encountered in the analysis, let us begin by studying the model 
restricted to boxes Qi C A. We take the side-length £ of order 
£ ~ ej^^ without being completely precise yet. This scale is the 
fundamental length scale in the problem as revealed below in 
formula (|5]l. 

In a fixed box Qi, the first calculation we do is an optimiza- 
tion of the Hamiltonian —T-Lq^ assuming the spin variables have 
small deviations from a fixed direction. Here —Hq/, is the ana- 
log of ([T]i in Qd with free boundary conditions. 

For convenience, denote a^ '■— O-x — \Qi\^^ ^zeo ^^ ^^^ 



Eg, (a) := 



E 



A' 



Expanding —Hq,, to second order in the 9 variables. 



sup 
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(4) 

The first term has the typical order of magnitude e^f^ for d > 3 
and e^£^ log£ for d = 2 while the second term is typically of 
order ei'^/'\ 

Neglecting the second term, directions of presumed ordering 
are obtained by neglecting the second term and optimizing the 
RHS in ip. For ip fixed, the optimal choice for the deviation 
variables Ox is 6^ = ecos{ij:)gQ^^x where gq^^x = -A"^ • a^- 

In evaluating the validity of this calculation one has to worry 
about three issues. First, we must understand in what sense the 
assumption of small deviation from a fixed angle is valid. If we 
denote the Dirichlet energy of spin configurations in boxes Q^ 
by 



x-a|=l 

the only a priori control of small amplitude deviations we 
are able to obtain regarding Eq^ is that it is costly energeti- 
cally for 8Qf,{a) > 4e^ \Qe\- One cannot ask for more than 
this because approximate ground states, which take the form 
ax = {cos{egQ^^x), sm{egQ^^x)), have this Dirichlet energy. 

Second, the calculation itself imposes an upper bound on 
£: It is a fact that gq^.x has typical order of magnitude in di- 
mensions d = 2, 3, 4 of i, Vt, Vlog^ respectively. This indi- 
cates that naive expansions breakdown beyond the length scales 
e~^, e~^, exp (e~^), respectively since eg is of constant order 
If we take fluctuations into account, in two dimensions this im- 
poses £ < e^^. If one believes in ordering, nonlinear effects 
must play a role at larger scales. 

Conversely, and of particular significance in two dimen- 
sions, moderate deviations of the field may make the term 



OUT. 



zeQt 



relevant. In two dimensions, this imposing 

a lower bound on £ > e ^ ^ | log e | ^ ^ . There is j ust enough room 
between the two constraints to make sense of the computation. 



If we take £ = e ^ | log e| 2 t for 7 e (0, 1/2) fixed, we are 
able to keep \gx\ < | loge]"''' (even taking into account fluctu- 
ations) while O ( e X^zgQ^ Q^z ) is lower order with probability 
exponential in | log ej^"^'''. 

A Key Computation. Next we present a computation of both 
physical and technical interest which allows us to turn the naive 
analysis presented above into a rigorous Peierls argument. The 
idea is that low energy spin configurations are composed of 
a fast oscillating approximate ground state (like gqi^x above) 
and a long wavelength contribution. Subtracting off the fast os- 
cillating component renormalizes the Hamiltonian into a more 
tractable form. 

For any connected set R C Z'^ (not necessarily a box), 
we introduce the Dirichlet Laplacian A|^ and represent 
spin configurations on R in angular coordinates by (Jx = 
{cos{9x),au-i{9x)). Using the field g'^ = [-Ag + £-^]-^ax, 
we make the change of variables 4>x ^ 9x ~ £C0s{9x)g'x- This 
change of variables should be compared with the optimizer for 
fixed angle ijj described after (|4). The mass^ term £^^ is intro- 
duced in order to keep the g'^ small in modulus, as long as the 
field ax behaves in a typical way inside the region R. 

The Hamiltonian ([T]) restricted to R transforms according to 
-Hfl,(a|f7O)K^/C(0, 100) with 



/C((/),|0O) = ^[cos( 
(xy) 



i]4E 



mijcos^f 



(5) 



where rn^ = J2v\x-v\=i^9'v ~ 9x]'^- There are errors in mak- 
ing this transformation, but these errors can be controlled if 
Dirichlet energy of cr in i? is smaller than 4e^|i?| and the field 
a behaves in a typical way inside R. Note that this transfor- 
mation introduces the length scale e^^ in an explicit way - rrix 
is typically of order e^. Thus, phase boundaries separating the 
two pure phases 9^0 and 9 ^ n have width of order e^ ^ . 
The Modified Peierls Argument. Next we explain how to use 
(|5]l to implement a Peierls argument. Taking inspiration from 
II28II . we define a second scale 




so that £< e^^ < L. 

Given a spin configuration ct on A, a box Q^ C F is bad 
for (J if one of two things occurs. THe first possibility is that 
the Dirichlet energy in Qi is substantially larger than the en- 
ergy scale Eqj (a) set by the randomness in Qi (this is typi- 
cally 0{e^£'^). The second possible source of bad behavior, if 
the Dirichlet energy is smaller than this scale, is that the aver- 
age angle tpQi, associated with a in Qi is bounded away from 
{0, tt} by a chosen cutoff ^. 

A region of space F C A is called contour for a if it is 
maximally connected union of boxes {Ql{z) : z G LZ''} so 
that for each Ql, there is a bad cube Qi within distance 3L/2. 
The goal, as with any Peierls argument, is to show that large 
contours are unlikely to occur at low temperature. 

Given a spin configuration a and an associated contour F, to 
extract energy cost from the existence of F we compare a with 
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a new spin configuration a which agrees with either a or the 
reflection of a across the 62 axis on each component of AAr\r. 
We require that a is within some (5 ^ ^ of either or tt on the 
whole of r. 

We construct a in a few steps. First, we find a layer C sur- 
rounding r which itself has "thickness" of order L. C has two 
properties: at the boundary of C, \iJx • ei| > 1/2 for all x and 
also ax • ei is of constant sign on each connected component 
of C. In what follows, a contour will be called a ± contour 
depending on the sign of cr^: • ei on the component of C which 
separates F from 00. 

It is worth remarking that it is here that our restriction on 
the dimension enters. By definition, at the boundary of a con- 
tour, the Dirichlet energy Eq,, is under control as is the spin 
average, i^'^ X^^eQ '-'x^ ^'^^ ^^^ boxes Qp C C. The lower the 
dimension, the more strongly this control restricts the size of 
"defects". 

Having found the boundary layer £, we change to the vari- 
ables inside this boundary layer, fixing the spins in AAr\£. This 
will cost us energetically, but less than we ultimately gain by 
extracting energy from the bad behavior on the contour 

Next, we replace cr in £ by the spin configuration a' associ- 
ated with the optimizer of (|5]l subject to the boundary condition 
given by a on Kfq\C. In each connected component of C, one 
can show that, for typical realizations of disorder, a' is uni- 
formly close to one of ±ei as long as dist(a::, C^) > L/4. This 
follows from the form of dD as long as L ^ e^ . 

The final step is to modify a' as follows. We first optimize 
ICa{4>) in the domian 

A:={xeCLir : dist(x, AAr\(£ U F)) > L/2} 

with free boundary conditions. Typically the optimizer 0_4 is 
close to throughout A. We want to "glue" a' in A'^ to the 
spin configuration 77 determined by (J>_a with minimal cost in 
energy. To do this, on each internal component of AAr\yl we 
replace, as necessary, a' by its reflection a" across the 62 axis 
so that a'^ ■ ei has the same sign on all components of C\A. a 
is defined as a" on A7v\^ and is 77 (or its reflection) on A. 

The big hurdle in all of these considerations is to make sure 
regions where fluctuations of the disorder invalidate the above 
reasoning are sparse. This is rather delicate in two dimensions. 
Details may be found in ifTol . 

Conclusion. In this letter we addressed a 25 year old contro- 
versy regarding randomness induced ordering, showing rigor- 
ously that the random field can actual select ordered phases in 
two and three dimensions. The results demonstrate that the low 
dimensional qualitative behavior of such models agrees with 
the first work on the subject ||2] and elucidates contradictory 
results in the papers II12II26I . 

There are two ways to understand these results in two dimen- 
sions. On the one hand, if one thinks along the lines of Ander- 
son localization, the randomness localizes the spin wave and 
vortex dipole excitations which enable the Kosterlitz-Thouless 
phase. On the other hand, if one thinks about the homoge- 
nization theory of classical fields, the behavior we discuss is a 
consequence of imposing the hard constraint that spins lie on 
the unit sphere. 



Our work gives a very precise picture for the behavior of 
these systems at low temperature, bounds on transition temper- 
atures in terms of the strength of the randomness and tools to 
extend the stated results to a broad class of classical and quan- 
tum systems. 
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